Filaggrin and cytokines in respiratory samples of preterm infants at risk for respiratory viral infection

Respiratory viral infections (RVIs) are frequent in preterm infants possibly inducing long-term impact on respiratory morbidity. Immune response and respiratory barriers are key defense elements against viral insults in premature infants admitted to Neonatal Intensive Care Units (NICUs). Our main goals were to describe the local immune response in respiratory secretions of preterm infants with RVIs during NICU admission and to evaluate the expression and synthesis of lung barrier regulators, both in respiratory samples and in vitro models. Samples from preterm infants that went on to develop RVIs had lower filaggrin gene and protein levels at a cellular level were compared to never-infected neonates (controls). Filaggrin, MIP-1α/CCL3 and MCP-1 levels were higher in pre-infection supernatants compared to controls. Filaggrin, HIF-1α, VEGF, RANTES/CCL5, IL-17A, IL-1β, MIP-1α and MIP-1β/CCL5 levels were higher during and after infection. ROC curve and logistic regression analysis shows that these molecules could be used as infection risk biomarkers. Small airway epithelial cells stimulated by poly:IC presented reduced filaggrin gene expression and increased levels in supernatant. We conclude that filaggrin gene and protein dysregulation is a risk factor of RVI in newborns admitted at the NICU.

www.nature.com/scientificreports/ Relevant clinical data are displayed in Table 1, showing no differences between the RVI (n = 26 infants) and the negative control (n = 23 infants) group regarding the most significant neonatal clinical outcomes (P > 0.05 for all comparisons). Nine out of the 26 infants (35%) who suffered RVI associated clinical signs of respiratory stagnation or deterioration, as defined in the methods section.
The type and frequency of respiratory viruses detected were: hRV (n = 7; 27%), hCoV (n = 5; 19%), hAdV (n = 6; 23%; one in coinfection with hBoV and another with picornavirus hPIV), hBoV (n = 3; 12%), RSV (n = 4; 15%; one of them with hMPV) and one case of enterovirus. We further classified patients regarding the virus type, but due to the small sample size of each virus and data from each group, we cannot observe any statistical difference for any of the clinical parameters (Table 1).
Cytokine profile of NPA supernatants is different in infants who will not develop a RVI at NICU compared who those who are infected. When comparing the samples of virus-negative controls against negative pre-infection infants, we observed that several molecules were reduced (Fig. 2), including filaggrin (19.3 ± 21.5 vs. 29.8 ± 27.3 ng/mL; P = 0.02), MIP-1α/CCL3 (1.6 ± 2.0 vs. 7.7 ± 19.7 pg/mL; P = 0.03), MIP-1β/ CCL4 (1.2 ± 6.0 vs. 10.9 ± 26.3 pg/mL; P = 0.04) and MCP-1 (23.1 ± 102.7 vs. 136.7 ± 335.0 pg/mL; P = 0.02). When we classified subjects of study accordingly to sex, there is no clear pattern of cytokine differentiation, showing no sex-dependent differences ( Supplementary Fig. S1). Moreover, we classified samples from viral positive subjects at the time of infection, and we could not observe any difference related to virus difference for any of the molecules analyzed ( Supplementary Fig. S2), being sample size and data from each group a limitation.
We went on to analyze whether the differentially expressed cytokines could be used as markers for predisposition to viral infection, or as markers for viral infection itself. As can be observed in Table 2, the levels of filaggrin, MIP-1α/CCL3 and MCP-1 could be used as moderate predisposition biomarkers with AUC values > 0.68, with the following values of specificity (> 84%) and sensitivity (> 53%). In addition, the levels of filaggrin, VEGF, MIP-1β, RANTES/CCL5, HIF-1α and IL-1β were identified as infection biomarker molecules with AUC values > 0.72, having the next values of specificity (> 79%) and sensitivity (> 55%) ( Table 2).
Moreover, we performed logistic regression multivariable models where we combined the biomarker value of the different measured cytokines, and we observed that the combination of the categorical values of filaggrin, MCP-1 and IL-1β (applying Youden index cutoff values) is a good biomarker for viral infection predisposition, having an AUC of 0.72, with 68% of sensitivity and 85% of specificity, that can help in differentiating neonates with RVI risk at NICU (Table 2). Similarly, the logistic regression model using continuous values of filaggrin, RANTES/CCL5, MIP-1β/CCL4 and MCP1 has an AUC of 0.80, with 74% of sensitivity and 90% of specificity for discrimination of those neonates who are infected at the moment of analysis at the NICU ( Table 2). The graphical representation of the ROC curves for both individual molecules and logistic regression models are provided in the Supplementary Fig. S3.
The immune response of infants with RVI changes with the course and evolution of the infection. Different  www.nature.com/scientificreports/ the RVI infants (P > 0.05), which could reflect a permanent dysregulation of this protein in RVI infants compared to controls. Finally, HIF-1α was increased in negative post-infection samples compared to negative pre-infection samples (542 ± 1697.5 vs. 145.4 ± 584.0 pg/mL; P = 0.04; Fig. 2), while the opposite trend was observed for MCP-1 (41.3 ± 111.4 vs. 136.7 ± 335 pg/mL; P = 0.02). A non-significant decrease in MIP-1α/CCL3, MIP-1β, RANTES/ CCL5 and MCP-1 (Fig. 2) was observed in negative post-infection samples, which could be related to the resolution of the inflammatory response.
Filaggrin cellular gene expression frequency and protein levels are altered in the NPA of neonates who will in time became infected by RVIs in NICU. To obtain information about the basal state (at the moment of NICU admission) of the airway barrier, we compared negative pre-infection samples from RVI and from virus-negative control groups. Having observed in our previous results that free filaggrin levels are different in samples depending on if the neonate becomes infected or not, we next studied the frequency of neonates who basally (never infected and pre-infection) do express filaggrin gene (FLG) at NPA cellular level, in order to get an insight of its expression in the cellular compartment of the airways. Moreover, we analyzed filaggrin protein quantity in the cells obtained from the NPAs using western Blot. We also studied the frequency of infants expressing other relevant genes such as TLR7, CAPN14 and DSG1 to get more insight into the state of the barrier. As depicted in Fig. 3a, the frequency of neonates with airway expression of FLG was reduced in the pre-infection samples from the RVI positive (55%) compared to the negative control (74%) group; the opposite occurred for TLR7, where higher rates of expression were found in the RVI positive pre-infection (71%) than in the negative control (37%) group. Neither of these trends reached statistical significance (P > 0.05), although these differences may provide information about how is the expression of immune and barrier genes in the free , and negative post-infection samples from virus-positive premature neonates (n = 19) measured by ELISA or Luminex technology. Inter-group analysis was performed using Kruskal-Wallis test with multiple comparisons followed by Dunn's uncorrected post-test. Bar graphs represent median and interquartile ranges. *P < 0.05; **P < 0.01 and ***P < 0.001. www.nature.com/scientificreports/

Epithelial cells from healthy donors increase immune-gene expression and release higher filaggrin after viral stimuli plus wounds.
To study lung barrier status we used a RVI in vitro model of SAEC, using poly:IC (20 ng/mL) as viral stimuli 24 . We observed a reduction in FLG gene expression in SAEC with (relative expression = 0.19 ± 0.19) or without (relative expression = 0.07 ± 0.26) wounds compared to the untreated control; (P = 0.004 and P = 0.02 respectively; Fig. 4a). Regarding the cell culture supernatant, the opposite effect was observed, with increased, although non-significantly, levels in secreted filaggrin both in poly:IC treated (31,333 ± 37,915 ng/mL; P > 0.05) and in poly:IC plus wound (28,772 ± 30,216 ng/mL; P > 0.05) compared to controls (5950 ± 15,047 ng/mL) and with wounds (14,528 ± 20,039.4 ng/mL) (Fig. 4B).

Discussion
This study, carried out in a population of preterm infants admitted to NICU, shows a distinct pattern of immune profile and epithelial barrier gene expression in baseline NPA from neonates who go on to develop RVI. Our findings indicate that vulnerability to infection may be related to immunological response or epithelial barrier integrity. We would like to highlight that the control group were randomized samples, obtained over the same time as the RVI positive newborns in order to obtain more representative data from this group (Fig. 1). In spite of this, we believe that modulation of the control group's molecular and immune response to be unlikely due to the absence of RVIs. On analysis of respiratory samples from never infected infants (controls) compared to pre-infected neonates (preterm infants who went on to develop RVI) we observed several differences in immune response cytokines and other molecules (including filaggrin, MIP-1 [α/β] and MCP-1), as well as in the expression and location of filaggrin in NPAs, that may indicate a barrier-immune deconfiguration predisposing certain individuals to future infection. Differences between samples from controls and neonates with ongoing , AREG (n = 7-9), CAPN14 and DSG1 (n = 6-9) (relative expression as 2 −∆∆Ct ) measured by qPCR and (b) filaggrin (pg/mL) protein synthesis measured by ELISA in SAEC treated with 20 ng/mL of poly:IC (polyinosinic:polycytidylic acid), wounded, or treated with poly:IC and wounded for 24 h (n = 5-9). Inter-group analysis was performed using the Kruskal-Wallis test with multiple comparisons test followed by Dunn's post-test. Bar graphs represent median and interquartile ranges. SAEC small airway epithelial cells. *P < 0.05; **P < 0.01. www.nature.com/scientificreports/ viral infection in several molecules (filaggrin, VEGF, RANTES (CCL5), HIF-1α, IL-17A and IL-1 β) were also observed, indicating an immune response associated with viral infection. By means of in vitro experiments we replicated reduction of FLG expression in airway epithelial cells after exposure to viral particles, observing that the filaggrin released into the extracellular space may be related to a failure in a correct epithelization, leading to a vicious circle of barrier disruption and infection.
The lung's immune response begins at the respiratory epithelium. Although our results are not statistically significant, we observed a modification of the frequency of gene expression and quantity of filaggrin levels in cells from the NPA of premature neonates who went on to suffer respiratory viral infection. Filaggrin is a protein that maintains the epithelial barrier of the nose and skin, and filaggrin defects have previously been associated with spontaneous atopic dermatitis, acting as disease biomarkers, and with secondary lung inflammation in mice 25,26 . Filaggrin loss of function mutations have also been associated with airway diseases, with higher allergen sensitization in children 27 and early wheezing 28 .
Interestingly, we found differences in the expression of several immune-barrier molecules (filaggrin, MIP-1α/ CCL3, MIP-1β/CCL4 and MCP-1) in NPA supernatants of virus-negative controls compared to negative preinfection samples of RVI neonates. These differences point to a possible predisposition to infection in RVI infants characterized by modified levels of filaggrin, MIP-1α/CCL3, MIP-1β/CCL4 and MCP-1. Filaggrin, MIP-1α, and MCP-1 levels may be useful for predicting the risk of subsequent RVI during NICU admission, using the NPA obtained early after birth. As previously mentioned, filaggrin's role is to maintain the epithelial barrier, and its release into the extracellular space can be associated with the barrier's disruption, and the consequent increase in permeability to substances and microorganisms. MIP-1α/CCL3 is necessary for anti-viral immune responses, and as previously described is a vital molecule in the recruitment of the leukocytes as monocytes, T lymphocytes (CD8 + and CD4 + ) and B lymphocytes to the site of infection 29 together with MIP-1β, which is also involved in mucosal humoral immunity, in the trafficking of T lymphocytes through the vascular endothelium and in the production of cytokines performed by T helper cells, like the release of IL-6 and IL-4 30 ; while and MCP-1 has also been described as a key antiviral molecule, as seen in studies unraveling the importance of this molecule in the severe cases of COVID-19, where increase of this molecule was predictive of respiratory failure 31 .
Several immune compounds were identified in NPA supernatants of infants with ongoing RVIs which differed to those of virus-negative controls, including IL-1β, a cytokine that promotes antiviral immunity through the activation of the inflammasome complex with nuclear factor-κB, inducing immune cell trafficking, activation and cytokine production of TNF-α and IL-6, that are characteristic of antiviral defense 12,32 ; VEGF and HIF-1α, both related to BDP development, as HIF-1α is an upstream regulator of VEGF, and together with NF-κβ are involved in apoptosis, alveolar structure and lung capillary density in neonatal development, factors that may differentiate between healthy and RVI predisposes neonates 33 ; and finally, filaggrin, associated with epithelial barrier diseases, as previously described, being a marker of decrease of this molecule in the airway epithelium linked to higher permeability to viruses by barrier disruption 28 . Some of these molecules are indeed potential biomarkers for differentiating between infected and uninfected neonates in the NICU, as seen in the ROC curve analysis and in the logistic regression models where cytokines and molecules related to lung integrity are good biomarkers to differentiate subjects and to predict risk.
The onset of infection provokes a series of changes in the immune compound composition of the NPA. During and after infection, a further increase in VEGF and HIF-1α was observed, which may reflect enhanced angiogenesis, and possible are also related to the homeostasis of oxygen during lung infections 34 . Equally higher concentrations of IL-1β, IL-17A and RANTES/CCL5 were detected. IL-17A is involved in viral clearance in the airways, being secreted by T cells, which creates an inflammatory milleu that favors the action of neutrophils and NK cells 35 ; and RANTES (CCL5), that causes viral replication control by the antiviral characteristics of eosinophils and T cells 13,36 . During viral immune responses, it seems that a general increase of antiviral and proinflammatory molecules like MIP-1α/CCL3 and MIP-1β/CCL4 occurs, a process that is coherent as a mechanism for virus control and lung homeostasis at the first stages of infection 30 . Interestingly, MCP-1 was reduced in the ongoing and post-infection RVI group, as described after immune recovery in bronchiolitis 37 , as a mechanism that might be related to the finalization of inflammation and antiviral responses 10 .
Several studies have described that umbilical cord blood samples from preterm infants showed downregulation of immune pathways, including complement cascades in infants and upregulation in maternal blood, which may imply that besides lung barrier integrity, immune responses might be reduced in preterm neonates 38,39 . Regarding the molecular sensors of innate immunity in premature infants, previous publications describe that TLR4 and TLR2 receptors are reduced in animals, impairing neutrophil recruitment, while neutrophils from neonates themselves show impaired migration, phagocytosis and cytotoxic activity [40][41][42] . Adaptative immune responses may influence the risk of viral infections in newborns as their CD4 + T cells have delayed synthesis of IFNγ and IL-2 in presence of viral infections 40,43 . Genetic studies can be very helpful in the determination of risk factors, such as mannose-binding lectin (MBL) mutations causing MBL deficiency, which are associated with higher risk of sepsis in preterm neonates 44 , or increased IFI27 expression, associated with higher RSV infection severity 45 .
These findings from the literature suggest that the risk of RVI in premature neonates is a complex combination of alterations in epithelial barrier integrity and immune system development. Our study provides more insight into this subject, describing a profile of NPA-specific biomarkers that could be of use in prevention and risk management of newborn ICU respiratory viral infections.
In the in vitro model of RVI we observed that, just as in NPAs from premature neonates, viral stimulus provokes a reduction in filaggrin expression, and an increase of released levels, perhaps indicating that viruses can cause barrier disruption by preventing cornification and airway homeostasis 46  www.nature.com/scientificreports/ conditions 24 , and in infants with bronchiolitis where this Th2 environment can lead to asthma and allergic diseases in the long term 47 . Interestingly, epithelial barrier integrity is affected by reduced AREG expression, a protein previously associated with tissue integrity, thus having relationship with the increase in permeability of the airway epithelium 48 . Nevertheless, we could not imitate the predisposition of infants to viral infection by wounds, as disruption of SAEC monolayer did not affect the expression of the selected genes, indicating that, in future, other models for barrier disruption would be more appropriate. The limitations of this study consist, firstly, in the lack of airway epithelial cell samples from neonates before and after infection, to prove a direct link between the cellular models and the nasopharyngeal aspirate samples. We used a cellular model consisting in small airway epithelial cells, which are anatomically different from nasal epithelial cells, and therefore might not respond in the same manner to viral stimuli. We also randomized timing of control samples during NICU admission to obtain a group that is representative of all three states (preinfection, active infection, and post-infection), which might mask temporal differences in the control group. Another limitation resides in the small sample size of our study, mainly due to the difficulty of obtaining samples from premature infants. Nevertheless, we feel that our results are of sufficient importance to inaugurate filaggrin in neonates as a field of research.
In summary, this study shows that preterm infants who will go on to suffer respiratory viral infection during NICU admission have a distinct respiratory epithelium immune profile at their NPA supernatants, consisting of higher baseline concentrations of immune and barrier-related molecules which can be used as possible biomarkers for infection risk for preterm neonates at NICU. In addition, the dysregulation of filaggrin synthesis and release from the epithelium, found in NPA and in SAEC with wounds and Poly:IC highlight the possible role of this molecule in epithelial barrier disruption, favoring RVI and triggering chronic diseases in neonate infants admitted in the NICU.

Methods
Study population. The patients included in this study are part of a prospective observational respiratory viral infection surveillance study of preterm infants (below 32 weeks of gestation) from the Departments of Neonatology of La Paz University Hospital (Madrid, Spain) 4 and the Severo Ochoa University Hospital (Leganés, Madrid, Spain) between December 2018 and December 2020 and whose primary aim was to determine the viral etiology of respiratory infections during the NICU admission, the outcome of infected infants and the risk factors related.
The present substudy is part of the previously mentioned prospective study and has the aim to describe the local immune response and the expression and synthesis of lung barrier regulators in respiratory secretions of neonates who develop RVIs during the NICU; and in vitro models. Inclusion criteria were being a preterm neonate of less than 32 weeks of gestation who was included in the NICU. Patients older than 3 days at enrolment, those with major malformations, and those who died within the first week of life, were excluded. A flow chart with the recruitment and classification of the subjects of study is represented in the Supplementary Fig. S5. The study was approved by La Paz University Hospital Clinical Research Ethics Committee and Severo Ochoa University Hospital Clinical Research Ethics Committee (code number HULP-PI-2255). Informed consent was obtained from parents of the infants included in the study. All infants' data were treated anonymously and codified.
Sampling. Nasopharyngeal aspirates (NPA) for RVI surveillance and immunological analysis were collected within the first 3 days after birth and then weekly until discharge. Additional NPA were obtained in the event of respiratory symptoms or clinical suspicion of sepsis. RVI samples were selected by identification of infants with viral infection confirmed by PCR. For each participant, a first negative pre-infection NPA sample (taken when the infant is admitted at the NICU), a second positive sample taken during ongoing viral infection, and a third negative post-infection sample (obtained at least 1 or 2 weeks after the last positive NPA) were selected for the study. The control group was formed by randomized samples (we chose the samples from the non infected controls using samples from different times at each control NICU stay, to homogenize as much as possible) from premature infants admitted to the NICU with a similar gestational age to their RVI counterparts who did not develop RVI during hospitalization, collected at three different times (equivalent to the three samples taken from patients in the RVI group, and therefore happening in the same seasons). The selection of control samples at equivalent time points was carried out to reduce the possible bias caused by differences in gestational age between groups. The study design is represented in Fig. 1.
Clinical parameters were recollected for assessment of the neonate's state, including gestational age, birth weight, bronchopulmonary dysplasia, supplementary oxygen, oxygen at discharge, days to reach full enteral nutrition, length of stay; and finally clinical signs of respiratory stagnation or deterioration defined as no progress made in the withdrawal of respiratory support despite the passage of time, and without any other condition that justifies it and as increased need for respiratory support compared to the previous situation with no other cause, respectively.

NPA cytokine and chemokine assays.
Small airway epithelial cell culture. The expression and synthesis of lung barrier regulators in the study and control group were investigated in NPA. An in vitro model was developed to further examine the results obtained in vivo, using primary small airway epithelial cells (SAEC) from healthy subjects (Lonza, Basel, Switzerland) cultured in a specific growth medium (Promocell, Heidelberg, Germany). Media were supplemented with 100 U/mL penicillin and 100 µg/mL streptomycin and maintained at 37 °C in an atmosphere containing 5% CO 2 . Cells were incubated with 20 ng/mL of poly:IC (Sigma Aldrich, MO, USA). Disruption of epithelial barrier integrity was performed using a pipette tip.
Statistical analysis. Descriptive data were expressed as median and interquartile ranges for continuous variables, and percentages for categorical variables. Normal continuous variables were compared using ANOVA with the Bonferroni correction, or through T-tests; and the Mann-Whitney U test or the Kruskal-Wallis test with uncorrected Dunn's post-test was used for non-normal samples. Categorical variables were compared using a chi-squared test or Fisher's exact test. P-values < 0.05 were considered statistically significant. The Youden index was used to set the optimal cut-off point. Logistic regression models were performed to determine the possibility of using a combination of molecules in order to differentiate between conditions, using as variables either the continuous data, or the categorical classification using as cutoff values the ones determined by the Youden index. The Akaike information criterion (AIC) was used for the determination of the best fitting classification models. Analyses were carried out using Graph-Pad Prism 8 (San Diego, CA, USA), and multivariate logistic regression models were performed using R (https:// www.r-proje ct. org/).

Ethics statement.
Subject's parents signed the required informed consent. The study was conducted in accordance with the principles of the Declaration of Helsinki and was approved by the hospital's ethics committee.

Data availability
The data that support the findings of this study are available from the corresponding author, V.P., upon reasonable request. www.nature.com/scientificreports/